Report for Model 10

— for both jointly and separately

2024-08-10

1. Joint Modeling

1.1 Trace plots for convergence check

The current MCMC setting is:

  • 100,000 iteration;
  • 90,000 burn-in;
  • 10 thinning.

1.2 Gelman and Rubin’s convergence check

## Potential scale reduction factors:
## 
##         Point est. Upper C.I.
## HDevsum       1.01       1.04
## LDevsum       1.04       1.13
## dh0           2.13       4.31
## dh1           2.85       4.86
## dh2           1.81       3.05
## dl0           1.45       2.10
## dl1           2.50       4.11
## dl2           1.17       1.47
## dl3           1.82       3.04
## 
## Multivariate psrf
## 
## 2.54

1.3 ACF Plots

Here we plotted ACF plots for the following variables:

  • Total deviance;
  • Variables that didn’t pass the convergence check.

1.4 WAIC results

LevelH LevelL
DIC 1160.44074 22017.4891
DIC3 1123.74866 22087.2419
PWAIC 26.83821 187.8045
WAIC 1137.43899 22099.3759

2. Separate Modeling

2.1 Trace plots for convergence check

The current MCMC setting is:

  • 100,000 iteration;
  • 90,000 burn-in;
  • 10 thinning.

2.2 Gelman and Rubin’s convergence check

## Potential scale reduction factors:
## 
##         Point est. Upper C.I.
## HDevsum       1.01       1.03
## LDevsum       1.01       1.02
## dh0           7.61      14.54
## dh1           8.89      19.10
## dh2           3.34       7.49
## dl0           1.17       1.47
## dl1           1.18       1.47
## dl2           1.03       1.07
## dl3           1.14       1.37
## 
## Multivariate psrf
## 
## 8.35

2.3 ACF Plots

Here we plotted ACF plots for the following variables:

  • Total deviance;
  • Variables that didn’t pass the convergence check.

2.4 WAIC results

LevelH LevelL
DIC 1393.35314 27391.3106
DIC3 1245.95662 27461.9357
PWAIC 90.50368 121.9975
WAIC 1315.16366 27465.8546